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Absorption of energy quanta generated by quantum point contact results in the inelastic current 
through the double quantum dot placed nearby. In contrast to a single quantum dot, the inelastic 
current through the double quantum dot is sensitive to the energy dependence of the quantum 
point contact transmission, which can explain the experimentally observed features. We calculate 
the inelastic current as a function of microscopic parameters of the circuit. 



Double quantum dot (DQD) has been recently pro- 
posed as a detector of nonequilibrium noise generated 
by nearby mesoscopic devices [1|. This idea was experi- 
mentally realized in measurements of the nonequilibrium 
noise spectrum of a quantum point contact (QPC) de- 
tected by DQD US]- The experiments provided a lot of 
interesting and in some respect puzzling results. 

The noise detection is based on the generation of in- 
elastic current through DQD assisted by absorption of 
energy quanta emitted by QPC. To implement the noise 
measurement, DQD is brought into the Coulomb block- 
ade regime with the highest energy electron localized in 
one of its dots hereafter referred to as dot 1. The ex- 
change of the electron between the two quantum dots 
brings DQD into the excited state, the excitation energy 
A being fixed by the gate voltages. Absorbing an energy 
quantum, an electron tunnels from the low energy state 
in quantum dot 1 to the excited state localized in quan- 
tum dot 2. The tunnel barrier between the quantum dots 
is tuned to be much higher than the barriers between the 
dots and the adjacent leads, so that after each interdot 
tunneling event the electron almost immediately escapes 
into the adjacent electron reservoir. Another electron 
occupies quantum dot 1, and the system returns to the 
ground state, the unit of charge having been transferred 
through DQD [2|. The generated current is therefore pro- 
portional to the noise power on the excitation frequency 
of DQD. The nonequilibrium noise is generated by QPC, 
which is brought in a strongly nonequilibrium transport 
regime by application of transport voltage. At the same 
time, the plunger voltage applied to QPC controls its 
transmission. 

Theoretical calculations of the generated inelastic cur- 
rent have been performed in Ref. where it has been 
related to the nonequilibrium noise power Si(A/h) gen- 
erated by QPC at the frequency A/H. This noise power 
is given by the local current fluctuations in an arbitrary 
spatial point of QPC, 
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Based on the energy conservation law one concludes that 
increasing the QPC transport voltage Vqpc, the current 



through DQD will start at the point Vq PC = A/e, when 
the quanta with energy A appear in the nonequilibrium 
noise spectrum [4| . A puzzling feature of the experimen- 
tal measurements is the independence of the threshold 
voltage Vq PC of the DQD excitation energy A for a fi- 
nite range of energies, contrary to the expectations based 
on the energy conservation law 0. 

In this Letter we provide a theoretical description of 
QPC-DQD system in the nonequilibrium regime that al- 
lows to relate experimental measurements of the inelastic 
current to microscopic parameters. We show that since 
DQD is an object extended over both sides of QPC, the 
local noise power (|T|) is not the relevant quantity for the 
inelastic DQD current. Rather, the noise power absorbed 
by DQD includes spatially nonlocal correlations of cur- 
rent fluctuations at positions of two quantum dots. The 
relevant voltage power is given by 



Sv(uj) = ( Z{x\, uj)SI(xi, uj) + Z(X2,U})SI(X2,U) 



(2) 

where Z(xi,uj) is the spatially dependent trans- 
impedance of the circuit, Xi denotes the position of the 
corresponding quantum dot, and u> = A/h is the absorp- 
tion frequency. In the case of the symmetric coupling of 
both quantum dots to QPC, which is called a symmetric 
circuit in what follows, the trans-impedance becomes in- 
dependent of spatial coordinate. Then the inelastic cur- 
rent is given by the expression similar to that obtained 
in Ref. ]l| 
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but with the nonlocal current noise power Si{lj) following 
from 

2 „oo 

S!(lu)=J2 dTeiWT {si(xi,T)Si(xj,0)) . (4) 
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The presence of nonlocal current correlations substan- 
tially modifies the resulting noise spectrum. So, while 
the energy dependence of the QPC transmission is not 
essential for the local current fluctuations Q]) [l], 2, [fl, 
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it becomes crucial for the spatially nonlocal fluctuations 
of the current ([4]). As the result, the generated inelas- 
tic current becomes sensitive to the energy dependence 
of QPC transmission. In particular, for the symmetric 
circuit the inelastic current turns to zero when the QPC 
transmission is energy independent within the transport 
voltage window. This can explain the independence of 
the DQD current onset voltage Vq PC on the DQD exci- 
tation energy A seen in experiment [2j. These findings 
are illustrated in Fig. [TJ If the QPC transmission am- 
plitude exhibits a plateau in its energy dependence, the 
current onset voltage is determined by the width of the 
plateau as long as that width exceeds the DQD excitation 
energy A (solid lines). In contrast, the threshold voltage 
equals A/e for the continuously rising QPC transmission 
(dashed lines). 

To emphasize the importance of the energy dependence 
of QPC transmission for the power absorbed by DQD, it 
is advantageous to define an analogy of the Fano factor 
for a finite frequency noise .F(w) = Si(uj)/(2elQpc)- This 
value relates the power generated at frequency u> and 
defined by JU) to the direct current through QPC. Its 
explicit expression reads 
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where t e and r e are the transmission and reflection am- 
plitudes of QPC at energy e. The Fermi distributions 
in the left (source)/right (drain) reservoirs of QPC are 
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denoted as f e ■ Their chemical potentials differ by the 
QPC transport voltage and can be written as Hl/r = 
±eVQPc/2, the chemical potential in the unbiased QPC 
being taken as zero. For zero temperature, and for the 
absorption frequency w = A/ft Eq. ([5|) simplifies to 



F(A/h) = e(eV-A)- 
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It is evident from (J6j) , that the Fano factor is determined 
by the energy dependence of QPC transmission. The de- 
pendence of the finite frequency Fano factor on the QPC 
transport voltage is shown in the inset to Fig. [TJ The 
energy threshold of the Fano factor corresponds to the 
width of the plateau in QPC transmission. Furthermore, 
the Fano factor drops if the applied voltage is larger than 
the energy interval for the onset of the conducting chan- 
nel. In that case, the DQD inelastic current reaches sat- 
uration while the direct QPC current continues to grow. 
Using the Fano factor ([5|), the expression for the gener- 
ated inelastic current can be cast in the form 

Ibqd = 2eI QPC T(A/h) r DQD , (7) 
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FIG. 1: (Color online) Current through DQD as a function 
of QPC bias voltage. Solid lines correspond to a plateau-like 
energy dependence of the QPC transmission amplitude i e , 
dashed lines correspond to a continuous t e (see upper inset). 
For A smaller then the width of the plateau Eo in t £ , the 
threshold voltage is defined not by A but by Eo (the current 
onset for the first two solid curves almost coincides). For 
continuous t e the current onset is determined by eVcjpc = A 
(see the dashed curves). Parameters of the interaction ll{ = 
Ul — 1.0. Other parameters are relevant for the experiment 
in Ref. |2|]. Lower inset: dependence of the finite frequency 
Fano factor on Vqpc- 



where Tdqd = ^hrq 1S ^ c DQD rectification factor, to 
denotes the tunnel amplitude between the two quantum 
dots, and Rq = h/e 2 is the quantum resistance. 

In what follows we introduce the theoretical model for 
the coupled QPC - DQD system and outline the deriva- 
tion of the presented results. Since the maximal inelastic 
current is observed when a new conducting channel is 
opening in QPC, we concentrate on a single conducting 
channel of QPC. We distinguish two species of electrons 
in QPC, namely those coming from the right and the 
left reservoirs, and we describe them by the fermion field 
operators ?Pr/l( x ) 0- The electrons of each sort are in 
equilibrium with its own reservoir. Taking the position 
of the QPC potential barrier at x = 0, we represent the 
field on each side of it as 

dk 



^(x) = j |_ ^j) L (k) e l(pF+k)x + r e e- i{pF+k)x 
+ip R (k)t* e e- i ( pF+ V x \ , for x < 0, (8) 

+ip L (k)t £ e l(pF+k)x } , for x > 0. (9) 



Here pf denotes the Fermi wave vector at zero transport 
voltage. 

The structure of the interaction between DQD and 
QPC-channel plays a crucial role for the generation 
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of inelastic current, determining the trans-impedance 
Z(xi,uj). Due to the presence of external electrodes, the 
effective interaction becomes screened and time-retarded. 
Moreover, the presence of QPC violates the spatial ho- 
mogeneity of the interaction. Therefore, we can write the 
interaction term in the action in the form 

Aim = -e 2 dx dtdt'Ui(\x\, \x - Xi\,t- t') x 



i=l,2 



h(x, t)fii(t') 



(10) 



Here h\ and fi2 are the particle number operators in each 
quantum dot, and h(x) is the operator of density fluctu- 
ations in the conducting channel at point x. We assume 
that the quantum dots are situated far away from the 
potential barrier of QPC, one on each side of it (see ex- 
perimental setup of Ref. Q). In the detection regime, 
the total occupation of DQD is fixed to ni(t) +ri2(t) = 1. 
This allows us to use a pseudo-spin 1/2 description of 
DQD. We associate the states localized in the quantum 
dots 1 and 2 with the spin-up and spin-down states re- 
spectively. The charge transfer between the two quan- 
tum dots corresponds to the spin-flip between the ground 
state spin-up and the excited state spin-down. The inter- 
action term (|10p can be splitted into the interaction with 
the total charge of DQD, and the interaction with the z- 
component of DQD pseudo-spin, S z = ni—ni- Note that 
the product tt — eS z , is proportional to the operator of 
the DQD dipolc moment. Omitting the interaction with 
the total charge of DQD, the relevant interaction can be 
represented as the one between the dipole moments of 
DQD and QPC, which in the Fourier-transformed form 
reads 
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where the QPC dipole moment P{lo) is defined as 
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The forward and backward inelastic scattering ampli- 
tudes are the only relevant ones for the one-dimensional 
motion of electrons in QPC. They are given in terms of 
Fourier transforms J dxe~ ig ( x ~ Xi 'Ui(\x\, \x — Xi\) at wave 
vectors \q\ <^p F and q = ±2p F . We assume that the in- 
teraction is strongly screened, and it takes place only in 
a small region of the size of the screening length around 
each quantum dot. Then the behavior of scattering am- 
plitudes at small wave vectors is smooth, and we can 
approximate Ui(\q\ <C p F ) ~ Ui{q = 0) = U[ for the for- 
ward scattering. For the backward scattering we obtain 
U%(±2p F ) = U l (±2p F )e^ 2tpFX ' Taking into account the 
finite size of a quantum dot, one has to integrate over 
Xi within that size, which greatly diminishes the back- 
scattering amplitude because of the rapidly oscillating 



factors e ±2ippXi . On that account we neglect the back- 
scattering amplitudes in what follows. 

Using Eqs. dSJ) , d9j) , (fT2l . we represent the total dipole 
moment operator P in terms of the fields of right- and 
left-reservoirs as follows 



/V! = - V / ^ w x^+ {k ')UkW-e-u). 

(13) 

The combined effect of the scattering by QPC potential 
barrier and interactions with quantum dots is captured 
by the effective inelastic scattering amplitudes * be- 
tween the two species of fermion fields. Further calcu- 
lation shows that only the inelastic scattering between 
different species contributes to the noncquilibrium noise 
power. The amplitude w LR is given by 
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and the amplitude w RL is obtained from w LR by ex- 
change U[ <-> — {//• Using the representation of the 
electric current operator in the basis (JSJ, ([9]) we can 
express the dipole moment P, in terms of the Fourier 
transform of the current operator at frequency to 0, 0] , 
P(u>) = J2i=i 2 Z(xi 7 uj)I(xi,uj) introducing the spatially 
dependent trans-impedance Z(xi,u>) = 



Iff. At this 



point it becomes evident that the dipole moment in- 
teracting with DQD involves spatially nonlocal correla- 
tions of QPC current. In the case of symmetric circuit, 
ll{ = [//, the trans- impedance becomes independent of 
coordinate. Its expression in terms of the elements of 
effective electric circuit is provided in Ref. []]]. 

The generated inelastic current is calculated perturba- 
tively in the lowest order of QPC-DQD interaction em- 
ploying the Keldysh technique 13, U | • The total action 
is given by A — ^o + ^int, with the interaction part given 
by (fTTj) and the free part 



(15) 



Ao = f dt J ^x G x^x + $D " 1$ \ 



In the free part of the action, we used the semifermionic 
representation of the pseudospin degrees of freedom of 
DQD [HEJ]. Here $ = 0_) T is the spinor Grass- 
mann field. Each spin component <\> a is in turn a two- 
component field in the retarded-advanced space. D is 
the semifermionic Green's function, its retarded compo- 
nents are given by D R a (e) = [e — a A/2 + io]^ 1 , D R _ a = 
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D R + ~ D R _). The semifermionic spin-1/2 represen- 



tation imposes special rules of calculus in Keldysh for- 
malism. The Keldysh component is parameterized as 



D K {e) = D R (e)F s (e) 
F s (e) = tanh(^) ± 

pansion, each diagram is calculated taking F s once with 



F s {e)D A (e) with the function 
jtx- In the diagrammatic ex- 
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FIG. 2: Diagrammatic representation: a) Vertex of QPC dipol 
moment P; b) Vertex of DQD dipol moment -if; c) DQD cur- 
rent vertex 7dqd; d) Diagram for DQD current in the second 
order of interactions. II denotes the QPC polarization op- 
erator, and F is the DQD rectification factor. Dashed lines 
denote the semifermionic Green's functions D(e); Solid lines 
denote the Green's functions in QPC. 



the plus and once with the minus sign, and the half-sum 
of the results is taken at the end of the calculation . 

The Green's function of the fermions from the reservoir 
X in QPC channel is given by G^, L (e, k) = [e =p vpk — 
■to] , the Keldysh component is obtained by usual 



rules 



The energy e is counted from the Fermi 



level of the unbiased QPC. 

The operator of the current through the double quan- 
tum dot can be written in the pseudospin representa- 
tion as /dqd = —etoS y /h. In the Keldysh formalism, its 
quantum component is proportional to the matrix a x act- 
ing in the retarded-advanced space /qq D = o x ® /dqd ■ 
The current is calculated using a diagrammatic repre- 
sentation derived from the Keldysh path integral with 
the action (Tl5|) (see Fig. [2]). The diagram for the cur- 
rent in the second order of perturbation in the interac- 
tion is shown in Fig. [2ji), the corresponding analyti- 
cal expression reads /dqd = Tr|ff[|. The polariza- 
tion operator fl = (P(u)P(— lu)^ can be represented 
in terms of current operators and trans-impedances as 
LI = /| J2i=i 2 Z(x%, u))I(xi, w)| 2 ^. For the symmetric 

circuit we obtain TL(oj) = \Z{lo)\ 2 Si{uj) with the noise 
power Si(uj) given by (|1]). Explicit calculation of the di- 
agram shown in Fig. [5Ji) leads to the final expression for 
DQD current ©. 

The suppression of the DQD current for Vqpc exceed- 
ing the energy conservation threshold A/e represents a 
profound feature of the current voltage characteristics of 
the double dot system under consideration, which has 
no analogy in the shot noise induced current through a 
single quantum dot [141 ]. It can be understood in simple 
terms, if one considers the generation of DQD current as 
a kind of a Coulomb drag experiment. Indeed, the di- 
agrammatic representation for the inelastic current Fig. 
[2jl) is almost identical to the diagramms for the drag 
current, the difference being the presence of a nonequi- 



librium polarizatio n op erator II for QPC instead of the 
current operator j{| [l5| ■ To realize the drag, the particle- 
hole (p-h) symmetry has to be violated in both compo- 
nents of the system. Its violation in the gated DQD is 
explicit. In the case of symmetric circuit, the violation of 
the p-h symmetry in the noncquilibrium QPC is possible 
only if the transmission amplitudes at energies differing 
by A are different. Indeed, the initial bosonic excitation 
that provides the energy to the DQD is an electron-hole 
pair in QPC with energy A. The Hamiltonian describing 
the propagation of such a pair through QPC is p-h asym- 
metric only for different transmission amplitudes for the 
electrons and the holes. 

The author is grateful to L. Glazman, A. Kamenev, 
and A. Levchenko for illuminating discussions. The 
author appreciates and enjoyed the hospitality of the 
William I. Fine Theoretical Physics Institute, Univer- 
sity of Minnesota, where a part of this work has been 
performed. The author acknowledges Financial sup- 
port from DFG through Sonderforschungsbercich 508 
and Sonderforschungsbercich 668. 



[1] 
[2] 
[3] 



[7] 

[8] 

[9] 
[10] 

[11] 



[12] 

[13] 
[14] 



[15] 



R. Aguado and L. P. Kouwenhoven, Phys. Rev. Lett. 84, 
1986 (2000). 

V.S. Khrapai, S. Ludwig, J. P. Kotthaus, H.P. Tranitz, 
W. Wegscheider, Phys. Rev. Lett. 97 176803 (2006). 
S. Gustavsson, M. Studer, R. Leturcq, T. Ihn, K. En- 
sslin, D. C. Driscoll, A. C. Gossard, Phys. Rev. Lett. 99, 
206804 (2007) 

Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

G. B. Lesovik, JETP Lett. 49, 683 (1989). 

Note that index R/L relates to the reservoir. It should 
not be mixed with chirality that is not conserved in the 
presence of potential scattering. 

The one-dimensional dispersion relation e = ftv^ |fe| al- 
lows to replace k by e. 

Thierry Martin, les Houches Session LXXXI, edited by 

H. Bouchiat et. al. (Elsevier 2005) 

A. Levchenko, A. Kamenev larXiv:0809.1670l (2008'). 

L. V. Keldysh, Zh. Eksp. Theor. Fiz. 47, 1515 (1964); 

[Sov. Phys. JETP 20, 1018 (1965)]. 

A. Kamenev, Many-body theory of non- equilibrium 
systems, in Nanophysics: Coherence and Transport, 
H. Bouchiat, et al. (editors); pp. 177-246, Elsevier, Am- 
sterdam, 2005. 

V. N. Popov and S. A. Fedotov, Zh. Eksp. Teor. Fiz. 94, 
183 (1988) [Sov. Phys. JETP 67, 535 (1988)]. 
M.N.Kiselev, Int. J. of Mod. Phys. B, 20, 381 (2006). 
E. Onac, F. Balestro, L. H. Willems van Beveren, U. 
Hartmann, Y. V. Nazarov, and L. P. Kouwenhoven, 
Phys. Rev. Lett. 96, 176601 (2006). 

L. Zheng and A. H. MacDonald, Phys. Rev. B 48, 8203 
(1993); Alex Kamenev and Yuval Oreg, Phys. Rev. B 52, 
7516 (1995); K. Flensberg, B. Y.-K. Hu, A.-P. Jauho, 
and J. M. Kinaret, Phys. Rev. B 52, 14761 (1995). 



